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Dissolvable  films  of  silk  fibroin  for  ultrathin 
conformal  bio-integrated  electronics 

Dae-Hyeong  Kim  and  Jonathan  Viventi  ef  al.* 

Electronics  that  are  capable  of  intimate,  non-invasive  integration  with  the  soft,  curvilinear  surfaces  of  biological  tissues 
offer  important  opportunities  for  diagnosing  and  treating  disease  and  for  improving  brain/machine  interfaces.  This  article 
describes  a  material  strategy  for  a  type  of  bio-interfaced  system  that  relies  on  ultrathin  electronics  supported  by  bioresorbable 
substrates  of  silk  fibroin.  Mounting  such  devices  on  tissue  and  then  allowing  the  silk  to  dissolve  and  resorb  initiates  a 
spontaneous,  conformal  wrapping  process  driven  by  capillary  forces  at  the  biotic/abiotic  interface.  Specialized  mesh  designs 
and  ultrathin  forms  for  the  electronics  ensure  minimal  stresses  on  the  tissue  and  highly  conformal  coverage,  even  for  complex 
curvilinear  surfaces,  as  confirmed  by  experimental  and  theoretical  studies.  In  vivo ,  neural  mapping  experiments  on  feline  animal 
models  illustrate  one  mode  of  use  for  this  class  of  technology.  These  concepts  provide  new  capabilities  for  implantable  and 


surgical  devices. 

Strategies  for  bio -integrated  electronics  must  overcome 
challenges  associated  with  the  mismatch  between  the  hard, 
planar  surfaces  of  semiconductor  wafers  and  the  soft,  curvi¬ 
linear  tissues  of  biological  systems.  These  differences  in  mechanics 
and  form  lead,  almost  invariably,  to  low-fidelity  coupling  at  the 
biotic/abiotic  interface  and  limited  long-term  tissue  health.  The 
difficulties  are  most  pronounced,  and  the  solutions  are  perhaps 
most  important,  in  systems  designed  for  brain/computer  interfaces 
(BCIs).  State-of-the-art  penetrating  microelectrode  arrays  consist 
of  sharp  shanks,  typically  10  x  10  arrays  of  pins  with  base  widths 
~80pm,  lengths  ~1.5mm  and  pitch  ~400pm  (ref.  1).  These 
arrays  are  rigid  and  inflexible  because  of  their  construction 
from  blocks  of  silicon,  which  also  supports  their  conventional 
wafer-based  electronics.  They  are  valuable  for  research  in  BCIs, 
but  they  damage  tissue  and  do  not  offer  long-term  electrical 
interface  stability2  because  of  unwanted  biological  responses  to 
the  electrodes.  Comparable  BCI  performance  can  be  achieved 
with  non-penetrating,  surface  electrode  systems  that  are  minimally 
invasive  and  provide  greatly  improved  stability3-5  with  minimized 
inflammation.  Standard  clinical  subdural  electrode  arrays  are  useful 
for  BCIs  (ref.  6)  but  their  widely  spaced  (~lcm),  large  contact 
electrodes  (~0.35  cm  diameter)  spatially  undersample  the  electrical 
signals  present  on  the  surface  of  the  brain7.  Decreasing  the  spacing 
and  size  of  the  measurement  points  can  improve  BCI  performance 
by  providing  access  to  high  temporal  and  spatial  frequency  signals8. 
Such  designs,  however,  demand  excellent  conformal  coverage  over 
the  highly  convoluted  brain  surface,  to  ensure  direct  coupling 
between  brain  tissue  and  the  electrodes. 

Reducing  the  thickness  of  the  substrate  decreases  the  bending 
rigidity,  thereby  improving  conformal  contact.  Unfortunately, 
clinical  arrays  and  even  the  thinnest  devices  designed  for  research 
have  thicknesses  (700  pm  and  >  10  pm  (refs  9,10),  respectively)  that 
are  larger  than  desired  to  ensure  conformal  contact.  Analogous 
systems  based  on  stretchable  substrates  have  also  been  explored 
in  other  neural  interfaces11,12  but  typically  with  similar  or 
larger  thicknesses.  In  conventional  designs,  ultrathin  geometries 
(that  is,  <10  pm)  are  impractical,  because  the  films  are  not 
sufficiently  self-supporting  to  be  manipulated  effectively  during 
fabrication  or  implantation. 


Silk  is  an  appealing  biopolymer  as  a  temporary,  soluble 
supporting  substrate  for  this  application  because  it  is  optically 
transparent13,14,  mechanically  robust  and  flexible  in  thin-film 
form15-17,  compatible  with  aqueous  processing18,19  and  amenable  to 
chemical  and  biological  functionalization13,20.  The  silk,  in  both  the 
non-treated  and  methanol-treated  formats,  is  biocompatible21,22, 
bioresorbable23  and  water  soluble  with  programmable  rates  of 
dissolution15,16.  Moreover,  recent  work  demonstrates  the  ability 
of  silk  films  to  serve  as  a  platform  for  transistors23  and  various 
classes  of  photonic  devices24,25.  The  process  of  preparing  silk 
substrates  for  the  purposes  reported  here  began  with  material 
derived  from  Bombyx  mori  cocoons,  and  followed  published 
procedures18,19.  Briefly,  boiling  the  cocoons  in  a  0.02  M  aqueous 
solution  of  sodium  carbonate  for  60  min  removed  sericin,  a 
water-soluble  glycoprotein  that  binds  fibroin  filaments  in  the 
cocoon  but  can  induce  undesirable  immunological  responses21,26. 
An  aqueous  solution  of  lithium  bromide  at  60  °C  solubilized 
the  fibres  and  subsequent  dialysis  removed  the  lithium  bromide. 
Centrifugation  followed  by  micro  filtration  eliminated  particulates 
to  yield  solutions  of  8-10%  silk  fibroin  with  minimal  contaminants. 
Casting  a  small  amount  of  the  solution  on  a  flat  piece  of 
poly(dimethylsiloxane)  (PDMS)  followed  by  crystallization  in  air 
(~12h)  yielded  uniform  films  (thickness  of  20-50  pm)  (Fig.  la) 
that  were  subsequently  removed  from  the  PDMS  for  integration 
with  separately  fabricated  electronics. 

For  the  systems  described  in  the  following,  ultrathin,  spin- 
cast  films  of  polyimide  (PI)  served  as  a  support  for  arrays  of 
electrodes  designed  for  passive  neural  recording.  Control  devices 
consisted  of  otherwise  similar  layouts,  but  formed  using  standard 
photolithographic  procedures  applied  directly  on  commercial  PI 
sheets  with  thicknesses  of  25  and  75  pm  (Supplementary  Fig.  SI). 
Anisotropic  conductive  film  (ACF)  bonded  to  electrode  pads  at 
one  end  of  the  arrays  provided  electrical  connection  to  external 
data  acquisition  systems  (Supplementary  Fig.  S2).  Ultrathin  PI 
films,  with  or  without  mesh  layouts,  cannot  be  manipulated 
effectively  for  processing,  interconnecting  or  implanting  onto 
the  brain  because  of  their  extreme  flexibility  and  mechanical 
fragility.  For  these  cases,  the  fabrication  process  exploited  layers 
of  PI  spin-cast  onto  silicon  wafers  coated  with  sacrificial  films  of 


*  A  full  list  of  authors  and  their  affiliations  appears  at  the  end  of  the  paper. 
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Figure  1 1  Schematic  illustration  and  images  corresponding  to  steps  for  fabricating  conformal  silk-supported  PI  electrode  arrays,  a,  Casting  and  drying 
of  silk  fibroin  solution  on  a  temporary  substrate  of  PDMS;  5-15-pm-thick  silk  film  after  drying  for  12  h  at  room  temperature,  b,  Steps  for  fabricating  the 
electrode  arrays,  transfer  printing  them  onto  silk  and  connecting  to  ACF  cable,  c,  Schematic  illustration  of  clinical  use  of  a  representative  device  in  an 
ultrathin  mesh  geometry  with  a  dissolvable  silk  support. 


poly  (methylmethacrylate)  (PMMA)  (left  frame  of  Fig.  lb).  After 
electrode  fabrication,  the  mesh  structure  devices  underwent  further 
etching  to  remove  unwanted  parts  of  the  PI.  The  processing  was 
completed  by  dissolving  the  PMMA  layer  with  acetone,  transfer 
printing  the  entire  assembly  to  a  film  of  silk  and  connecting  the 
ACF,  yielding  easily  manipulated  bioresorbable  neural  recording 
systems.  See  schematic  illustrations  and  images  in  Fig.  lb.  In 
all  cases,  the  arrays  consisted  of  30  measurement  electrodes 
(Au,  150  nm)  in  a  6  x  5  configuration,  each  with  dimensions 
of  500  pm  x  500  pm  and  spaced  by  2  mm.  Interconnection  wires 
were  protected  by  a  thin  (~  1.2  pm)  overcoat  of  PI  to  prevent 
contact  with  the  tissue  or  surrounding  fluids.  Choosing  the 
thickness  of  the  PI  passivation  layer  to  match  that  of  the  PI 
substrate  locates  the  interconnects  at  the  neutral  mechanical  plane, 
thereby  minimizing  the  potential  for  bending-induced  mechanical 
fracture.  Details  of  the  fabrication  steps  appear  in  the  Methods 


section.  Electrode  arrays  were  implanted  by  placing  them  on 
the  exposed  brain  (after  craniotomy)  and  then  flushing  with 
saline  to  dissolve  the  silk.  This  procedure  induced  spontaneous, 
conformal  wrapping  of  the  device,  as  illustrated  schematically  for 
the  mesh  design  in  Fig.  lc.  After  the  measurements,  the  electrode 
array  can  be  easily  removed,  because  of  the  attachment  of  the 
electrodes  to  the  ACF. 

The  sequence  of  images  in  Fig.  2a  shows  the  dissolution  process 
for  a  representative  case  (7-pm-thick  PI  film,  connected  to  ACF 
on  a  silk  substrate  with  a  thickness  of  ~25  pm)  inserted  into 
warm  water  (~35°C).  As  the  silk  substrate  disappears,  the  total 
bending  stiffness,  El ,  diminishes  markedly  because  of  its  cubic 
dependence  on  thickness.  Computed  results  appear  in  Fig.  2b  and 
Supplementary  Fig.  S3c  for  PI  thicknesses  of  2.5  and  7  pm.  To 
highlight  the  benefits  of  reduced  thickness,  the  inset  shows  the 
ratio  of  El  for  these  two  cases.  Through  programmed  control 
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Figure  2  |  Time-dependent  changes  as  the  silk  substrate  dissolves,  a,  Dissolution  of  the  silk  through  submersion  in  warm  water,  b,  Total  bending  stiffness 
of  7  pm  and  2.5  pm  electrode  arrays  on  supporting  silk  films  as  a  function  of  the  thickness  of  the  supporting  silk  film.  Inset:  The  ratio  of  bending  stiffness 
between  7  pm  and  2.5  pm.  c,  Time-dependent  change  in  the  volume  of  a  silk  film  during  dissolution  (left  frame)  and  bending  stiffness  calculated  for  silk 
treated  in  70%  ethanol  for  5  s  for  two  different  array  thicknesses  (right  frame).  The  5  s  ethanol  treatment  increases  the  dissolution  time  from  minutes  to 
about  1  h. 


of  the  dissolution  rate  by  modifications  of  the  silk  protein 
secondary  structure15,16,  these  changes  in  El  can  be  designed 
to  occur  over  periods  of  time  ranging  from  seconds  to  years, 
depending  on  requirements.  Figure  2c  shows,  as  an  example,  the 
dissolution  rate  of  silk  film  slightly  treated  with  ethanol  (left 
frame)  and  the  computed  time  dependence  of  El  in  devices  that 
employ  more  thorough  ethanol  treatment  (right  frame).  The  error 
range  for  silk  thickness  measurement  is  ±7%.  See  Supplementary 
Information  for  detailed  conditions.  This  dissolution  time  can  be 
lengthened  even  more  by  extending  the  treatment  time  to  days 
or  weeks15;  the  corresponding  time  dependence  of  El  appears  in 
Supplementary  Fig.  S4. 

To  examine  the  ability  of  these  systems  to  conform  to  relevant 
surfaces,  we  carried  out  experiments  using  a  human  brain  model, 
following  the  basic  steps  shown  in  Fig.  lc.  Figure  3  provides  images 
for  various  cases  after  washing  with  saline,  including  relatively  thick 
control  devices  that  do  not  incorporate  silk.  Clearly,  the  extent 
of  conformal  coverage  increases  with  decreasing  thickness;  the 
mesh  design  provides  further  improvements,  as  shown  in  Fig.  3d, 
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Supplementary  Figs  S5  and  S6.  To  reveal  the  underlying  mechanics, 
we  carried  out  systematic  and  quantitative  studies  on  well-defined 
surfaces  that  capture  certain  basic  features  of  the  curvature  of  the 
brain.  The  first  set  of  experiments  explored  wrapping  the  devices 
on  isolated  and  overlapped  cylindrical  surfaces.  Figure  4a  shows  the 
simplest  case  of  a  device  with  bending  stiffness  £/,  thickness  h,  width 
b  and  length  21,  wrapped  on  a  cylinder  with  radius  R.  Analytical 
expressions  for  El  can  be  written  for  the  multilayer  structures  of 
Fig.  1  in  terms  of  materials  properties  and  geometries,  as  described 
in  the  Supplementary  Information.  For  the  wrapped  state  to  be 
energetically  favourable, 


where  y  is  the  adhesion  energy  per  unit  area.  The  bottom  frame 
of  Fig.  4a  compares  the  above  relation  with  a  series  of  experiments 
(Supplementary  Fig.  S7).  The  data,  the  error  range  of  which  is 
±5%,  are  consistent  with  y  ~  10  mj  m-2,  which  is  comparable  to 
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Figure  3  |  Neural  electrode  arrays  of  varying  thickness  on  simulated 
brain  models  to  illustrate  flexibility,  a,  Schematic  illustration  of  trends  in 
thickness  and  structure  that  improve  conformal  contact,  b,  Series  of 
pictures  illustrating  how  the  thickness  of  the  electrode  array  contributes  to 
conformal  contact  on  a  brain  model,  c,  Magnified  view  of  the  pictures  in  b. 
d,  Image  of  an  electrode  array  with  a  mesh  design  on  a  dissolvable  silk 
substrate.  The  arrows  indicate  struts  in  the  mesh  that  help  to  stabilize  the 
Au  interconnects  after  dissolution  of  the  silk.  The  inset  illustrates  the  high 
degree  of  conformal  contact  that  can  be  achieved  on  the  brain  model  once 
the  silk  substrate  has  been  dissolved. 

reported  values  for  wet  interfaces27.  Reducing  the  thickness  provides 
clear  benefits,  for  example,  wrapping  cylinders  using  only  capillary 
adhesion  forces  is  possible  for  R  ~  1  cm  when  h  <~  15  pm. 

A  pair  of  overlapped  cylinders  represents  a  simple  model  for 
a  gyrus  of  the  brain.  Figure  4b  shows  cylinders  with  radii  R,  a 
centre-to-centre  separation  of  2d  and  connected  by  a  smooth  arc  of 


radius  r0,  at  the  angular  position  00  =  sin  1  [d/(R  +  r0)].  The  contact 
angle  of  a  thin  film  with  one  cylinder,  6 ,  can  be  shown  to  be 


Rsin0  dR0cos0  /  y  \  /  d  d0cos0\ 
d  —  RsinO  +  (d  —  RsinO)2  \yc  /\  RsinO^~  Rsin20  ) 


where  yc  is  given  in  equation  (1).  The  solution  of  equation  (2)  takes 
the  form  0  —  0(d/R,  y / yc ).  For  y  <  yc,  the  energy  has  a  minimum 
at  6  —  0,  and  the  film  does  not  wrap  around  the  cylinders.  Partial 
wrapping  occurs  to  a  contact  angle  of  0  (that  is,  contact  for  angles 
between  0  and  0  <  00)  for  yc  <  y  <  y'c,  where  y'c  is  obtained  from 
equation  (2)  with  6  =  00  as  given  in  the  Supplementary  Information. 
For  y  >  y',  wrapping  is  complete  (that  is,  conformal  contact  for 
angles  between  0  and  60).  By  comparing  equation  (2)  with  the 
experiment  in  Supplementary  Fig.  S8,  the  extracted  adhesion  energy 
per  unit  area  is  y  =  lOmJm-2.  Results  appear  in  the  bottom 
frame  of  Fig.  4b,  where  the  parameters  correspond  roughly  to 
features  on  the  brain  model:  R  =  6.14  mm,  d  =  5.93  mm  and 
r0  =  1.72  mm.  The  error  range  of  the  data  is  ±5%.  (Experimental 
images  appear  in  Supplementary  Fig.  S8.)  By  substituting  0  with 
0O  in  equation  (2),  the  critical  thickness  for  conformal  contact  is 
obtained  as  h0  —  4.9  pm  for  the  present  system;  that  is,  devices 
thinner  than  ~4.9  pm  achieve  conformal  contact  on  this  surface. 
The  experimental  results  are  consistent  with  this  calculation. 

Cylindrical  surfaces  such  as  those  of  Fig.  4a, b  are  developable; 
the  brain  is  not.  As  a  model  non- developable  surface,  we  examined 
the  case  of  a  hemispherical  substrate.  Figure  4c  shows  results 
for  electrode  arrays  with  sheet  designs  at  thicknesses  of  7  and 
2.5  pm  and  with  an  open  mesh  layout28  at  2.5  pm,  each  on  a  glass 
hemisphere  with  a  radius  of  curvature  of  6.3  mm.  With  only  water 
capillarity  as  the  adhesion  force,  the  mesh  electrode  array  achieves 
excellent  conformal  contact.  The  sheets  show  comparatively  poor 
contact,  with  large  wrinkles,  even  for  the  thinnest  case  (that 
is,  2.5  pm).  Mechanical  analysis  of  a  simple  model  reveals  the 
underlying  physics.  The  left  frame  of  Fig.  4d  shows  the  mechanics 
model  for  the  sheet  design,  which  consists  of  a  circular  film  with 
radius  r  +  w  wrapped  onto  a  sphere  with  radius  R.  The  central 
green  part  denotes  a  PI  plate  of  radius  r,  tension  stiffness  (£/z)PI  and 
equi-biaxial  bending  stiffness  (EI)m.  The  yellow  ring  corresponds 
to  a  multilayer  structure  of  PI  and  Au,  of  width  w,  tension  stiffness 
(Eh) composite  and  equi-biaxial  bending  stiffness  (El) composite-  For  the 
film  to  wrap  around  the  sphere,  the  required  minimum  adhesion 
energy  per  unit  area  is  obtained  analytically  as 


(ei) n  (Eh) pi  r  ( 

R 2  r2  Jo  V 


R  •  x\2 
—  sin  —  xdx 
x  RJ 


2  W  (El)  composite  w(Eh)  composite 

7r2  +  r 


(3) 


A  mechanics  model  for  the  mesh  design  appears  in  the  right  frame 
of  Fig.  4d,  which  consists  of  only  a  circular  strip  of  a  corresponding 
multilayer  of  PI  and  Au.  In  this  case,  the  minimum  adhesion 
energy  per  unit  area  is 


For  the  case  that  w  <^r,  ]/csheet  in  equation  (3)  is  always  larger 
than  ycmesh  in  equation  (4);  that  is,  ycsheet  >  ycmesh.  The  inference  is 
that  the  open  mesh  design  requires  a  much  lower  adhesion  energy 
than  the  corresponding  sheet,  thereby  leading  to  greatly  improved 
ability  for  conformal  coverage,  as  shown  in  Supplementary  Fig. 
S9a.  Supplementary  Fig.  S9c  shows  critical  adhesion  energies  for 
films  with  thicknesses  up  to  80  pm.  For  a  thickness  of  2.5  pm  and 
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w/r  =  4,  ycsheet  =  29.1  mjm-2  for  the  sheet,  which  is  more  than 
12  times  larger  than  the  mesh  ycmesh  =  2.4  mj  m-2.  In  addition,  the 
mesh  design  involves  membrane  strains  that  are  smaller,  by  roughly 
a  factor  of  w/r,  compared  with  sheets  with  similar  thickness. 
For  the  experimental  mesh  systems,  this  ratio  is  of  the  order 
of  1/4.  As  a  result,  for  a  representative  critical  wrinkling  strain 
of  0.1%,  nearly  two  thirds  of  the  sheet  will  wrinkle.  Under  the 
same  conditions,  the  entire  mesh  gives  perfect,  conformal  contact. 
Finally,  the  normal  (peeling)  interfacial  stress  for  the  mesh  is  only 
1  /4  of  that  for  the  sheet  (Supplementary  Figs  S9b  and  S9d),  leading 
to  improved  adhesion  and  reduced  forces  applied  to  the  substrate. 
See  Supplementary  Information  for  details. 

In  vivo  neural  monitoring  experiments  on  a  feline  animal 
model  demonstrated  the  practical  implications  of  these  favourable 
mechanics.  The  tests  involved  an  anaesthetized  cat  with  its  head 
fixed  in  a  stereotaxic  apparatus  and  its  eyes  focused  on  a  monitor 
that  subtended  28°  x  22°  of  space.  An  initial  craniotomy  and 
durotomy  exposed  a  2  x  3  cm  region  of  cortex.  The  electrode  arrays 
covered  much  of  the  visual  cortex  as  shown  in  the  left  frames 
of  Fig.  5a-c.  Visual  stimuli  consisted  of  full-field  drifting  gratings 
presented  for  Is  at  2 Hz  with  a  spatial  frequency  of  0.5 cycles 
per  degree.  Gratings  were  presented  at  two  different  directions 
over  eight  different  orientations  (16  unique  stimuli).  However,  the 
responses  obtained  from  all  16  unique  stimuli  were  averaged  to 
obtain  the  largest  possible  signal-to-noise  ratio. 

Three  kinds  of  electrode  array  were  used  for  comparison: 
76-pm-  and  2.5-pm-thick  sheets  and  a  2.5  pm-thick  mesh.  The 
second  two  included  dissolvable  silk  supports.  The  left  images  of 
Fig.  5a-c  illustrate  the  progressively  improved  conformal  contact 
with  reduced  thickness  (that  is,  76  pm  to  2.5  pm,  in  Fig.  5a  and  b, 
respectively)  and  with  introduction  of  the  mesh  (that  is,  Fig.  5c). 
The  right  frames  of  Fig.  5a-c  demonstrate  the  effectiveness  of 
decreasing  the  electrode  thickness  and  the  mesh  structure  on 
physiological  measurements  of  brain  activity. 

In  particular,  these  frames  show  the  average  evoked  response 
measured  at  each  electrode,  each  plotted  in  a  spatial  arrangement 
that  corresponds  to  the  images  in  the  left  frames.  Prominent 
visually  evoked  potentials  are  observed,  particularly  a  strong  PI 00 
response.  The  PI 00  response  is  a  ‘positive’  evoked  response 
typically  occurring  at  100  ms  after  the  stimulation  onset29.  The  PI 00 
responses  shown  in  Fig.  5  are  plotted  positive  down,  by  convention. 
The  background  colour  of  each  plot  illustrates  a  quantitative 
measure  of  the  evoked  response  signal  quality.  This  measure  of 
signal  quality  was  calculated  by  dividing  the  root  mean  square  (rms) 
amplitude  of  each  average  electrode  response  in  the  200  ms  window 
immediately  after  the  presentation  of  the  visual  stimulus  by  the 
rms  amplitude  of  the  average  1.5  s  window  immediately  preceding 
the  stimulus  presentation.  The  colour  bar  at  the  bottom  of  Fig.  5c 
provides  the  numerical  scale  for  all  of  the  colours  used  in  Fig.  5a-c. 
This  measurement  serves  as  a  quantitative  metric  of  the  electrode 
performance,  because  the  uniform  nature  of  the  stimulation  is 
expected  to  evoke  similar  responses  across  the  entire  visual  cortex. 

In  each  case,  28  of  the  30  electrode  channels  were  recorded  and 
evaluated  for  evoked  potential  response,  as  coloured  in  green  to  red. 
Two  channels,  indicated  in  grey,  were  reserved  as  local  references, 
as  required  by  the  recording  apparatus,  and  were  not  evaluated.  The 
channels  with  high  and  low  rms  amplitude  ratios  are  coloured  green 
and  red,  respectively.  The  76  pm  (Fig.  5a)  electrode  array  exhibited 
the  lowest  performance  with  a  mean  rms  amplitude  ratio  of  all  28 
channels  of  3.6  d=  1.8.  This  was  due  to  poor  contact  at  many  of 
the  electrodes.  The  2.5  pm  array  (Fig.  5b)  showed  better  conformal 
contact  and  correspondingly  a  higher  mean  rms  amplitude 
ratio  of  5.2  =b  3.9.  However,  the  higher  standard  deviation  and 
correspondingly  wide  spectrum  of  the  red  and  green  channels  on 
the  array  indicate  that  although  many  electrodes  recorded  excellent 
signals,  approximately  half  of  the  electrodes  still  had  poor  contact 


2L  4 


Figure  4  |  Mechanical  modelling,  theoretical  predictions  and  measured 
properties,  a,  A  thin  film  wrapped  around  a  cylinder  of  radius  R.  The 
unwrapped  and  wrapped  states  appear  in  the  top  and  centre  frames, 
respectively.  The  bottom  frame  compares  the  mechanics  model  and 
experiments,  b,  A  thin  film  wrapped  around  two  overlapped  cylinders.  The 
top  and  centre  frames  show  the  unwrapped  and  wrapped  states, 
respectively.  The  bottom  frame  shows  a  comparison  between  the 
mechanics  model  and  experiments,  c,  Images  of  electrode  arrays  (76  pm 
sheet  in  left  top,  2.5  pm  sheet  in  right  top  and  2.5  pm  mesh  in  bottom 
panel)  wrapped  onto  a  glass  hemisphere,  d,  Mechanics  models  for  sheet 
(left  frame)  and  mesh  (right  frame)  designs. 


with  the  brain  and  recorded  weak  responses.  The  2.5  pm  mesh 
electrode  (Fig.  5c)  showed  the  best  performance,  with  nearly  all 
channels  in  good  contact  and  a  still  higher  mean  rms  amplitude 
ratio  of  5.7  =b  3.0.  The  lower  standard  deviation  of  the  2.5  pm  array 
illustrates  that  most  of  the  electrodes  recorded  good  responses. 

Figure  5d  shows  representative  single -channel  data  from  one  of 
the  2.5  pm  mesh  electrodes.  A  sleep  spindle  is  observed  with  good 
signal  amplitude  and  signal-to-noise  ratio.  This  collective  set  of 
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observations  is  consistent  with  the  systematic  mechanics  studies 
described  previously.  We  did  not  observe  any  evidence  of  immune 
response.  Histology  data  from  related  types  of  device  implanted 
under  the  skin  exhibited  no  inflammation  after  4  weeks,  as  shown 
in  Supplementary  Fig.  SI 0. 

Although  purely  passive  electrode  systems  serve  to  demonstrate 
the  advantages  and  underlying  aspects  of  these  systems,  the 
same  approaches  are  compatible  with  fully  active  electronics  and 
optoelectronics.  This  technology  allows  intimate  integration  of 
finely  spaced  electrode  systems  with  living  tissue,  enabling  the  kind 
of  reliable  biotic/ abiotic  interface  with  moving,  biological  structures 
that  will  be  required  for  chronically  implanted,  high- resolution 
medical  devices.  This  improved  electrode/tissue  interface  has  the 
potential  for  a  positive  impact  on  human  health  in  many  modes 
of  use. 

Methods 

Thick  electrode  array  (>25  jxM)  fabrication.  Commercial  PI  films  (Kapton, 
Dupont)  with  thicknesses  of  25  and  75  pm  were  attached  to  a  temporary  carrier 
substrate  consisting  of  a  glass  slide  coated  with  PDMS.  After  cleaning  the  surfaces 
with  acetone,  isopropyl  alcohol  and  deionized  water,  electron  beam  evaporation 
formed  uniform  coatings  of  metal  (Cr/Au,  50/1450  A).  Photolithography  and 
patterned  etching  yielded  arrays  of  interconnect  lines.  Thin  layers  of  PI  (thickness 
~1.2  pm)  spin-cast  and  patterned  by  reactive  ion  etching  left  only  the  ends  of  the 
lines  exposed.  Further  deposition  and  patterning  defined  square  metal  electrode 
pads  at  these  locations.  Peeling  these  away  from  the  PDMS-coated  glass  slide  and 
bonding  them  to  an  ACF  cable,  using  procedures  described  in  a  separate  section, 
completed  the  fabrication.  Supplementary  Fig.  SI  provides  a  schematic  diagram 
and  images  of  the  process. 

Thin  electrode  array  (<10  pM)  fabrication.  The  fabrication  in  this  case  used 
a  carrier  silicon  wafer  coated  with  a  thin  (~1.2  pm)  spin-cast  layer  of  PMMA 
(A2,  MicroChem).  The  device  substrate  consisted  of  a  film  of  PI  (Sigma  Aldrich) 
spin-cast  onto  the  PMMA.  Procedures  similar  to  those  described  for  thick  devices 
formed  the  metal  electrodes  and  PI  overcoat.  After  fabrication,  the  ultrathin  devices 
were  released  by  dissolving  the  sacrificial  PMMA  layer.  Transfer  printing  with  a 
PDMS  stamp  delivered  the  devices  to  dry  silk  film  substrates,  coated  with  ~9%  silk 
solution  as  an  adhesive.  The  final  step  involved  bonding  of  an  ACF  cable. 

Mesh  electrode  array  (<10  pM)  fabrication.  The  first  and  last  parts  of  the 
fabrication  sequence  were  identical  to  the  steps  outlined  in  the  previous  section. 
The  only  difference  was  the  addition  of  a  step  to  remove  certain  regions  of 
the  polymer  layers  (that  is,  PI  and  underlying  PMMA)  by  oxygen  reactive 
ion  etching  through  a  mask  (design  in  Supplementary  Fig.  S3)  to  define  the 
mesh  structure.  Detailed  dimensions  are  as  follows:  thickness  ~2.5  pm,  contact 
electrode  size  500  pm  x  500  pm,  mesh  width  ~250  pm.  (See  more  details  in 
Supplementary  Fig.  S5.)  This  etching  immediately  followed  the  formation 
of  the  electrode  pads. 


;L 

200  ms 


ACF  connection.  The  contact  pads  on  the  electrode  array  were  first  aligned  with 
the  ACF  cable.  Metal  clips  were  used  to  apply  pressure,  spread  evenly  over  the 
contact  pad  area  using  a  piece  of  PDMS  inserted  between  the  ACF  and  the  clips. 
Next,  the  clamped  sample  and  ACF  were  placed  in  an  oven  preheated  to  ~150  °C 
for  ~15  min.  This  process  formed  a  strong  mechanical  bond  between  the  electrode 
array  and  the  ACF  with  low  electrical  resistance. 


Figure  5  |  Photographs  and  data  from  animal  validation  experiments. 

a-c,  Image  of  an  electrode  array  on  a  feline  brain  (left)  and  the  average 
evoked  response  from  each  electrode  (right)  with  the  colour  showing  the 
ratio  of  the  rms  amplitude  of  each  average  electrode  response  in  the 
200  ms  window  (plotted)  immediately  after  the  presentation  of  the  visual 
stimulus  to  the  rms  amplitude  of  the  average  1.5  s  window  (not  shown) 
immediately  preceding  the  stimulus  presentation  for  a  76  pm  (a),  2.5  pm 
(b)  and  2.5  pm  mesh  (c)  electrode  array.  The  stimulus  presentation  occurs 
at  the  left  edge  of  the  plotted  window.  In  all  three  images,  the  occipital  pole 
is  at  the  bottom  of  the  frame  and  the  medial  is  at  the  right.  The  scale  bars 
at  the  bottom  of  c  indicate  the  spatial  scale  for  the  left  frames  and  the 
voltage  and  timescales  for  the  right  frames  of  a-c.  The  colour  bar  at  the 
bottom  of  c  provides  the  scale  used  in  the  right  frames  of  a-c  to 
indicate  the  rms  amplitude  ratios,  d,  Representative  voltage  data  from  a 
single  electrode  in  a  2.5  pm  mesh  electrode  array  showing  a 
sleep  spindle. 


Data  acquisition  and  processing.  The  electrode  arrays  were  connected  to  a 
Neuralynx  DigitalLynx  data  acquisition  system  through  the  ACF  and  a  custom 
electrode  interface  board.  The  board  appears  in  Supplementary  Fig.  SI  1. 
Supplementary  Fig.  S2  shows  the  connected  electrode  array,  ACF  ribbon  and 
circuit  board.  The  signals  were  sampled  at  the  standard  DigitalLynx  sampling 
rate  of  32,556  Hz  per  channel.  The  high-pass  filter  was  set  at  1  Hz  and  low-pass 
filter  at  9  kHz.  Custom  MATLAB  software  (The  MathWorks)  was  used  for  offline 
processing.  Signals  were  down-sampled  to  2,713  Hz  and  further  low-pass  filtered  at 
50  Hz.  Responses  were  averaged  for  each  stimulus  and  plotted  per  electrode. 

Animal  experiments.  Experiments  were  conducted  in  accordance  with  the  ethical 
guidelines  of  the  National  Institutes  of  Health  and  with  the  approval  of  the 
Institutional  Animal  Care  and  Use  Committee  of  the  University  of  Pennsylvania. 
Surgical  and  stimulation  methods  were  as  described  in  detail  previously29-31 .  Briefly, 
adult  cats  (2. 5-3. 5  kg)  were  anaesthetized  with  intravenous  thiopental  with  a 
continuous  infusion  (3-10  mg  kg-1  h_1)  and  paralysed  with  gallamine  triethiodide 
(Flaxedil).  Heart  rate,  blood  pressure,  end-tidal  C02  and  electroencephalograph 
were  monitored  throughout  the  experiment  to  assure  depth  and  stability  of 
anaesthesia  and  rectal  temperature  was  kept  at  37-38  °C  with  a  heating  pad.  The 
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surface  of  the  visual  cortex  was  exposed  with  a  craniotomy  centred  at  Horsley 
Clarke  posterior  4.0,  lateral  2.0. 

For  visual  stimulation,  the  corneas  were  protected  with  contact  lenses  after 
dilating  the  pupils  with  1%  ophthalmic  atropine  and  retracting  the  nictitating 
membranes  with  phenylephrine  (Neosynephrine).  Spectacle  lenses  were  chosen 
by  the  tapetal  reflection  technique  to  optimize  the  focus  of  stimuli  on  the  retina. 
The  position  of  the  monitor  was  adjusted  with  an  x-y  stage  so  that  the  area 
centralae  were  centred  on  the  screen.  Stimuli  were  presented  on  an  Image  Systems 
(Minnetonka)  model  M09LV  monochrome  monitor  operating  at  125  frames  s-1  at 
a  spatial  resolution  of  1 , 024  x  786  pixels  and  a  mean  luminance  of  47  cd  m-2. 
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Supplementary  Information 

Silk  dissolution  test 

By  altering  the  secondary  structure  of  a  silk  film,  one  can  program  the  amount  of  time 
it  takes  for  the  film  to  dissolve  in  water.  For  this  paper  we  want  the  film  to  dissolve 
within  minutes,  or  within  hours.  No  treatment  was  required  to  make  the  film 
dissolve  in  minutes.  Creating  some  beta  sheet  structure  by  exposing  to  70%  ethanol 
for  approximately  5  seconds  increased  the  dissolution  time  to  approximately  1  hour. 
To  determine  a  dissolution  rate,  we  made  51  square  inch  films,  exposed  them  to 
ethanol  for  5  seconds,  put  them  in  a  room  temperature  water  bath  and  measured  their 
dry  weight  after  a  certain  time  in  the  water  bath.  See  Fig  2c. 


The  bending  stiffness  of  the  thin  film 

The  cross  sectional  geometry  of  the  thin  film  is  illustrated  in  Fig.  S3a.  There  are  n 
gold  bricks  (size  bmxhm,  Young’s  modulus  EAu  =78  GPa  and  Poisson’s  ratio 
v Au  =  0.44  )  surrounded  by  PI  (size  b  x  h  ,  Young’s  modulus  EPI  =  2.5  GPa  and 
Poisson’s  ratio  vPI  =  0.34 ).  The  distance  between  the  neutral  axis  and  bottom  of  the 
thin  film  is 
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where  h'  is  the  distance  between  bottoms  of  gold  bricks  and  thin  film.  The 


bending  stiffness  of  the  thin  film  is 
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The  bending  stiffness  of  the  thin  film  on  a  silk  backing  substrate 

The  cross  sectional  geometry  of  the  thin  film  on  a  silk  backing  substrate  is  illustrated 
in  Fig.  S3b.  The  silk  backing  substrate  has  a  thickness  H  and  Young’s  modulus 
Esuk  =2.8  GPa  .  The  distance  between  the  neutral  axis  and  bottom  of  the  thin  film  is 
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The  bending  stiffness  of  the  thin  film  is 

(\ 

-h2  -  hy0  +  yl 
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A  thin  film  wrapping  around  two  overlapped  cylinders 


The  beam  theory  gives  the  unwrapped  part  of  thin  film  (above  the  connecting  arc  in 
the  center  frame  of  Fig.  4b)  is  (part  of)  a  circle  of  radius  r  given  by 

d 


r  =  • 


sin# 


-R  . 


The  bending  energy  in  the  thin  film  is  obtained  as 
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The  adhesion  energy  is 
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The  total  energy  of  the  wrapped  state  is  the  summation  of  the  above  bending  energy 
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and  adhesion  energy, 
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where  yc  is  given  in  Eq.  (1). 

y'c  is  obtained  from  Eq.  (2)  with  0  =  6O  as  y'c  =  jl  +  (l  +  X)R2 -  A)r02  Jj  and 
A  =  K)d  R  +  r^(R  +  r„)  —  d  sin  [/?/(/? +  ro)^]j. 


Mechanics  model  for  mesh  design 

For  the  sheet  design  as  shown  in  the  left  frame  of  Fig.  4d,  the  thin  film  is  modeled  as 
a  plate.  The  central  green  part  denotes  the  PI  plate  of  radius  r,  tension  stiffness 
(Eh\=E„h  and  equi-biaxial  bending  stiffness  (El)pi  =  EPI(\  +  vPI)li  j\2 ,  where 
EP,  =  Epi  j ~  Vpi  j  is  the  plane-strain  modulus  of  PI.  The  yellow  ring  is  the 

sandwiched  composite  of  PI  and  Au  of  width  w,  tension  stiffness 
msite  =  Ep/h  +  [EAu  - Eri)hm  and  equi-biaxial  bending  stiffness 
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-  is  the  distance  between  the  neutral  axis  and 


bottom  of  the  thin  film,  and  EAu  =  EAuj [l  -  vAu  j  is  the  plane-strain  modulus  of  Au. 
The  total  energy  of  the  wrapped  state,  which  is  composed  of  the  bending  energy  and 
membrane  energy  in  the  thin  film  and  the  adhesion  energy  between  the  thin  film  and 
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sphere,  is  given  analytically  as 
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For  the  thin  film  to  wrap  around  the  sphere,  U sheet  ^  0 ,  which  gives  the  required 


minimum  adhesion  energy  per  unit  area  yscheet  in  Eq.  (3).  The  maximum 


circumferential  membrane  strain  is 
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The  maximum  interfacial  normal  (peeling)  stress  is  obtained  as 
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For  the  mesh  design  as  shown  in  the  right  frame  of  Fig.  4d,  the  total  energy  of  the 
wrapped  state,  which  is  also  composed  of  the  bending  energy  and  membrane  energy 
in  the  yellow  composite  ring  and  the  adhesion  energy  between  the  thin  film  and  the 


sphere,  is  given  analytically  as 
2  nrw 
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From  Umesh  ^  0  for  the  thin  film  to  wrap  around  the  sphere,  the  required  minimum 


adhesion  energy  per  unit  area  y™esh  in  Eq.  (4)  is  obtained.  The  maximum 


circumferential  membrane  strain  is 
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For  representative  geometrical  parameters  used  in  the  current  study,  e.g.  w  =  250  pm , 
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r=  1  mm  and  7?=  1 0  mm,  s2esh  =0.07%.  The  maximum  interfacial  normal  (peeling) 
stress  is  obtained  as 


(S14) 


Figure  Captions 

Figure  SI.  Electrode  array  fabrication  process  using  thick  PI  film  (Kapton, 
Dupont,  USA),  a,  attach  PI  film  to  PDMS  coated  glass,  b,  electrode  array  fabrication, 
c,  ACF  connection. 

Figure  S2.  Images  of  electrode  array  after  connection  of  ACF  and  circuit  board,  a, 

electrode  array  with  thin  (<  10  pm)  substrate  thickness,  b,  electrode  array  with  thick 
(>  10  pm)  substrate  thickness. 

Figure  S3.  Schematic  diagram  for  analytical  model  and  its  modeling  result,  a, 

Cross  section  of  the  neural  sensor,  with  geometrical  parameters  illustrated,  b,  Cross 
section  of  the  neural  sensor  on  a  silk  backing  substrate,  c,  The  bending  stiffness  of 
the  neural  sensor  of  thickness  7  pm  and  2.5  pm  on  a  silk  backing  substrate. 

Figure  S4.  Time  dependent  bending  stiffness  change  for  7  pm  and  2.5  pm 
electrode  array. 
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Figure  S5.  Design  of  mesh  type  electrode  array. 

Figure  S6.  Images  of  mesh  type  electrode  array  after  dissolution  of  silk  substrate. 

a,  on  glass  cylinder,  b,  on  human  brain  model. 

Figure  S7.  Wrapping  experiments  on  single  cylinder,  whose  diameter  is  a,  3.5cm,  b, 
1.3cm,  c,  0.4cm  and  d,  0.1cm. 

Figure  S8.  Wrapping  experiments  on  overlapped  cylinder,  a,  angled  view,  b,  side 
view. 

Figure  S9.  Modeling  results,  a,  The  critical  adhesion  energy  and  b,  the  normal 
(peeling)  stress  between  the  fdm  and  sphere  surface  for  sheet  and  mesh  designs,  c, 
The  critical  adhesion  energy  for  sheet  and  mesh  designs,  d,  The  normal  (peeling) 
stress  between  the  film  and  sphere  surface  for  sheet  and  mesh  designs. 

Figure  S10.  A  histological  section  of  a  ~100  micron  thick  film  embedded 
beneath  the  epidermis  of  a  mouse.  The  silk  implant  (B)  is  integrated  well  between 
the  epidermis/dermis  (A)  and  the  hypodermis  (C).  The  section  was  taken  4  weeks 
after  implantation  and  stained  with  hematoxylin  and  eosin.  The  absence  of  any 
significant  chronic  inflammation  indicates  that  the  silk  films  are  well  tolerated  over 
the  examined  time  period. 
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Figure  Sll.  Design  of  the  circuit  board  used  to  connect  the  ACF  ribbon  to  the 
Neuralynx  DigitalLynx  data  acquisition  system. 
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a 


Spin  coat  PDMS, 
attach  Kapton  film  on  glass 


b 


Device  fabrication, 
cut  edge  of  kapton 


c 


Detach  from  glass, 
connect  ACF  film 


Figure  SI 
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Figure  S2 
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Figure  S4 
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Figure  S5 
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Figure  S6 
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Figure  S7 
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Figure  S8 
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Figure  S9 
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Figure  S11 


NATURE  MATERIALS  |  www.nature.com/naturematerials 


19 


